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ABSTRACT 

We show that it appears possible for starburst galaxies, like the nearby NGC 
253, recently identified as a TeV source by the CANGAROO collaboration, to 
emit a significant amount of high-energy 7-rays and neutrinos through hadronic 
processes in their cores. We suggest that proton illumination of the inner winds 
of massive stars can be a viable mechanism for producing TeV 7-rays and neu- 
trinos without a strong MeV-GeV counterpart. The rich stellar content of the 
starbursts, with millions of early-type stars concentrated in the central regions, 
where collective effects of the stellar winds and supernovae can produce a signif- 
icant enhancement of the cosmic ray density, provides an adequate scenario for 
TeV 7-ray generation. Close starbursts are also found to be potential sources 
for km-scale neutrino telescopes, like ICECUBE, within reasonable integration 
times. 

Subject headings: galaxies: starburst, gamma rays: observations, gamma rays: 
theory, galaxies: individual (NGC 253) 



1. Introduction 

Very recently, the CANGAROO collaboration reported the detection of the nearby 
starburst galaxy NGC 253 at TeV 7-ray energies (Itoh et al. 2002). The purpose of this 
Letter is to show that it is possible for starburst galaxies, like NGC 253, to produce TeV 
photons and neutrinos through hadronic processes. 
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2. A plausible sceneirio 

Collective effects of strong stellar winds and supernova explosions in star-forming regions 
are expected to result in particle acceleration up to multi-TeV energies, or even beyond (e.g. 
Montmerle 1979; Cassc & Paul 1980; Bykov & Fleishman 1992a,b; Anchordoqui et al. 1999; 
Bykov 2001). Recently, the region co-spatial with the Cygnus 0B2 association has been 
detected as a TeV 7-ray source by the HEGRA telescope array (Aharonian et al. 2002). 
TeV cosmic rays (CRs) accelerated in the association might be responsible for the high- 
energy 7-ray emission through the hadronic illumination of some suitable target (Butt et al. 
2003). A nearby EGRET source (3EG J2033+4118) has also a likely stellar origin (Chen et 
al. 1996, Romero et al. 1999, Benagha et al. 2001). 

In a star burst, where the ambient density of CRs is enhanced by a local high rate of 
supernova explosions and the collective effects of strong stellar winds (e.g. Bykov 2001), 7- 
rays can be produced by interactions of relativistic protons with the rich interstellar medium. 
This interaction produces neutral pions that quickly decay into 7-rays. The latter can also 
be produced by leptonic processes, like inverse Compton interactions with the strong FIR 
field, and relativistic bremsstrahlung in the ambient gas. All these processes have been 
modelled with a set of reasonable parameters for NGC 253 by Paglione et al. (1996). The 
expected total 7-ray flux, however, is below the EGRET sensitivity hmit, in accordance with 
the non-detection of the galaxy in the MeV-GeV range (Bhattacharya et al. 1994, Blom et 
al. 1999). 

A source of 7-rays not previously considered in the mentioned analysis is the production 
at the base of the strong stellar winds of early-type stars. When the relativistic particles 
are accelerated in the stellar wind itself, for instance through multiple shocks produced by 
line-driven instabilities (e.g. White & Chen 1992), the 7-ray flux from pion decays is limited 
to MeV energies. However, if the star is immersed in an external bath of relativistic particles 
accelerated at larger scales, e.g. in the core of the starburst, TeV protons can penetrate the 
base of the wind to produce TeV 7-rays. The injection of MeV — GeV particles into the wind 
of the stars will be strongly attenuated by the modulation effects of the wind itself. 

The dense medium in which these 7-rays are produced is, nonetheless, transparent to 
7-ray propagation. The optical depth to pair production is r ~ nci?, where n is the photon 
number density, a is the photon-photon cross section, and R is the distance that the photon 
must travel to escape. In order to compute n, we need the photon energy density at the base 
of the stellar wind, U. We shall assume that target photons in the wind have typical energies 
of ~ 1 eV. Then, U = ?>L/{A'kR^ c) ~ 2 x 10^° eV cm~^, where we have taken, for order-of- 
magnitude estimates, the following typical values: ~ 5 x 10^^ cm and L ~ 10^^ erg s~^. 
Then, n — U/{l eV), and the optical depth is very small, r ~ 2 x 10~^. The computation of 
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T makes use of the 7 — 7 (1 TeV — leV) cross section, cr ~ 2 x lO"^*^ cm^ (Lang 1999). Once 
TeV photons are produced and escape from the wind region, they must yet leave the galaxy 
to reach the Earth. The starburst has a FIR field with U ~ 200 eV cm'^ (e.g. Paglione 
et al. 1996), and photons have to travel distances of about 70-100 pc to leave the active 
region of the galaxy. The optical depth is, in this case, r = (200 photons cm~^) x (2 x 10~^° 
cm^) X (100 pc) ~ 1.2 xlO~^, which is also small enough to secure that photons can escape 
from the galaxy. 



3. NGC 253 

NGC 253 has been described as an archetypal starburst galaxy by Ricke et al. (1980), 
and it has been extensively studied from radio to 7-rays (e.g. Beck et al. 1994, Paglione et al. 
1996, Ptak 1997). The supernova rate is estimated to be as high as 0.2 — 0.3 yr^\ comparable 
to the massive star formation rate, ~ O.IMq yr~^ (Ulvestad et al. 1999, Forbes et al. 1993). 
The central region of this starburst is packed with massive stars. Watson et al. (1996) have 
discovered four young globular clusters near the center of NGC 253; they alone can account 
for a mass well in excess of 1.5xlO^M0 (see also Keto et al. 1999). Assuming that the star 
formation rate has been continuous in the central region for the last 10^ yrs, and a Salpctcr 
IMF for 0.08-100 Mq, Watson et al. (1996) find that the bolometric luminosity of NGC 253 
is consistent with 1.5 xIO^Mq of young stars. Based on this evidence, it appears likely that 
there are at least tens of millions of young stars in the central region of the starburst. Many 
stars might be obscured by the large amount of molecular material, as in the case of the 
massive Galactic cluster Westerlund 1 (Clark and Negueruela 2002). 

Physical, morphological, and kinematic evidence for the existence of a galactic super- 
wind has been found for NGC 253 (e.g. McCarthy et al. 1987, Heckman et al. 1990). 
This superwind creates a cavity of hot (~ 10^ K) gas, with cooling times longer than the 
typical expansion time scales. As the cavity expands, a strong shock front is formed on the 
contact surface with the cool interstellar medium. The shock velocity can reach thousands 
of kilometers per second, and iron nuclei can be efficiently accelerated up to high energies 
~ 10^° eV by diffusive mechanisms operating at large scales. Indeed, NGC 253 has been 
proposed as the origin for some of the observed ultra-high energy CRs (Anchordoqui et al. 
1999, 2002). 

The CANGAROO collaboration did not give a precise photon spectrum for NGC 253. 

The reported data were only the integrated flux above 0.5 TeV (securing an lla confidence 
detection) and the possible range of power law indices F =(1.85, 3.75) (Itoh et al. 2002). 
We can summarize these results expressing the differential photon spectrum as: F^{Ej) ~ 



-4- 



B {E^/lTeV)~^ cm-2 s-^ TeV-^ where B = (3.8 ± 1.2) x lO'^^ for V = 1.8, or B ^ 
(3.2 ± 1.0) X 10~^^ for r = 3.75. It is worth noticing that the direct extrapolation of the 
spectrum at TeV energies down to the GeV band would yield an integrated value compatible 
with the EGRET upper limits for NGC 253 only for power law indices T < 2. Spectra softer 
than those with F = 2 at TeV energies should present a break at a few hundred GeV. 



4. Hadronic origin of the TeV emission from NGC 253? 

The differential 7-ray photon number distribution from neutral pion decays at the source 
is given by (e.g. Gaisser 1990) N^{E^) = Np{E^) ^4:71 app{E^) N 2Z^^'^^o/o:. Here, E^ is the 
photon energy, = H/irip is the column density, is the proton mass, and a is the proton 
spectral index, such that Np = KpE'"^. In addition, ~ 1.5 is a correction factor that takes 
into account possible effects of heavier nuclei, app{E^) ~ 35 mb is the hadronic interaction 
cross section, and Z^'^^o ~ 0.17 stands for the fraction of kinetic energy of the parent proton 
that is transferred to the neutral pion during the collision. For a photon/proton spectral 
index a ~ 2 we get N^{E^) = 1.12 x lO'^^KpE'^ cm'^ TeV-^ (notice that 1 erg ~ 1 TeV). 
Assuming typical grammages for the stellar winds of OB stars in the range E ~ 50 — 150 g 
cm-2 (e.g. White 1985), we get N^{E^) = (3.4 - lO.l)^^^;-^ cm'^ TeV"^ . The luminosity 
in the CANGAROO observing range is 

/~20TcV 
N^{E^)E^dE^ (1) 
'0.5TeV 

= (4.7 - 13.9) X lO^^R^Kp erg s-\ 

where R is the size of the base of the wind and c is the speed of light (all numerical coefficients 
are for cgs units). The CR energy density is, on the other hand, 

/~20TeV 
Np{Ep)EpdEp = 9.9Kp erg s"^ = ^c^cr,©, (2) 
'IGeV 

where q is the amplification factor of the CR energy density with respect to the local value, 
<^CR,0 ~ 10~^^ erg cm~^. This imphes that the TeV luminosity can be written as L = (0.47 — 

lA)R'^q erg s~^. The latter corresponds to the luminosity of a single star illuminated by TcV 
CR protons. In the central region of the starburst, millions of stars could be illuminated at 
the same time. The total luminosity would then be Ltotai = — 1-4) i?f erg s^^, 

where N^^ is the number of contributing stars, and we have allowed for the possibility that 
not all of them will have the same R, nor will they be located within regions having the 
same CR enhancement factor. 
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Fig. 1. — Left (a): Results of 1000 Monte Carlo trials to compute the total TeV luminosity 
of a starburst lying at the distance of NGC 253. Right (b): The hadronic model expectation 
of neutrinos from NGC 253 compared with the upper limit for the neutrino background in 
a 1° X l°-generation neutrino telescope. 

If we now assume the following ranges for the different parameters involved: 1 x 10^^ < 
R[cm\ < 1 X 10^^ (Lamers 1999), 300 < <; < 3000 (Suchkov et al. 1993), and 10^ < A^^ < 
5 X 10^ (see above), there is ample room for the total luminosity to reach values above 10^^ 
erg s~^, i.e. the observed CANGAROO luminosity above 0.5 TeV, given the distance to 
NGC 253 (2.5 Mpc, de Vacoulers 1978). To show this, we have carried out Monte Carlo 
numerical simulations where the number of stars illuminated by the accelerated protons, the 
grammage, the size of the base of the wind for each star, and the ambient enhancement of 
CRs are independently allowed to take values within their assumed ranges. For each of these 
cases, the total luminosity was computed. We show the results in Fig. la, for a sample of 
1000 trials. Only for less than 40% of the trials, the total luminosity is below 10^^ erg s~^. 
This scenario then appears to be a plausible explanation for the TeV observations of NGC 
253. 

Note that multi-TeV electrons in the starburst central region will experience inverse 
Compton losses with the FIR fields (typical energy density ~ 200 eV cm~^) in the Klein- 
Nishina regime, which result in a very steep 7-ray spectrum (that can be approximated 
by power-law indices > 3.2 (although notice that the actual spectrum is not a pure power 
law, see e.g. Georganopoulos et al. 2001). The total leptonic contribution from the central 
starburst region at TeV energies is then negligible in comparison to the hadronic one. 
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MeV — GeV particles cannot easily penetrate the stellar winds, and hence these winds 
will glow mainly at TeV energies. The modulation of CRs by stellar winds is a complex 
subject that has been studied in detail only in the context of the relatively weak solar wind. 
In such case, three-dimensional models including diffusion and the effects of the terminal 
shocks have been developed (e.g. Kota & Jokipii 1983, Jokipii et al. 1993). These models 
can explain fairly well the basic observed features of the nucleonic CR component in the solar 
system, where strong gradients in the CR density are known to exist from interplanetary 
spacecraft measurements (e.g. Quenby et al. 1990).^ In OB stars, with extremely supersonic 
winds and mass loss rates orders of magnitude higher than the solar values, these effects 
should be far more pronounced, leading to an almost complete rejection of externally injected 
MeV — GeV particles at the base of the winds and thus to an absence of a strong (above 
EGRET sensibility) 7-ray counterpart at these energies. The modulation effect depends on 
the parameter e ~ uR^/D, where u is the stellar wind velocity, i?„ is the radius of the wind 
shell, and D is the diffusion coefficient, e measures the ratio between the diffusive timescale 
of the particle to their convectivc timescale. The former should be shorter than the latter for 
a TeV particle to enter into the wind cavity. For a typical O star, R^, > 10 pc and u ~ 2500 
km/s. Hence, the diffusion coefficient of the TeV particles should be at least 10^^ cm^ s~^, 
similar to what is estimated for our Galaxy (e.g. Berezinskii et al. 1990). 



5. Neutrino signal 

If the 7-ray emission is explained by the decay of neutral pions, their charged counter- 
parts must produce a neutrino flux. Following Alvarez-Muniz & Halzen (2002), the z/^ + i^^ 
neutrino flux F^{E^) produced by the decay of charged pions in the source can be derived 
from the observed 7-ray flux F^{E^) by imposing energy conservation (see also Stecker 1979): 

/ ' E^F^{E^)dE^ I E,F,(E,)dE, . (3) 

Here, (^™^) is the minimum (maximum) energy of the photons (neutrinos) and the 

pre-factor K — 1.^ The maximum neutrino energy is fixed by the maximum energy of 



•^Although it is highly dependent of the solar activity and the inclination angle of the incident particles, 
the proton flux decreases to less than ~ 30% of its initial value at 80 A.U. for 2.5 GeV particles, and to less 
than ~ 10% for 250 MeV particles. 

"^Note that there are two muon-neutrinos (out of the pion-to-muon-to-electron decay chain) for every 
photon produced. Each of the latter carries, on average, an energy whereas each photon carries on 
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the accelerated protons which can in turn be obtained from the maximum ob- 

served 7-ray energy {E';^'''') as: E^'''' ^ 6E^^ , E^'''' ~ 1/12E^'^ (Alvarez-Muniz 
& Halzen 2002). The minimum neutrino energy is fixed by the threshold for pion pro- 
duction, E"™"^ = 5 ^'^"^p^"^^^ — where 5 is the Lorentz factor of the accelerator rela- 

' P 2mp ' 

five to the observer. Assuming a Lorentz factor of order 1, the //-spectrum results in 
F^{E^) ~ 4.2 X 10-^2(E/TeV)-^-^5 TeV"^ cm'^ s^^ for the harder photon index, and 
F^{E^) ~ 9.5 X 10-^''^(^/TeV)-3-75 XeV"^ cm'^ for the softer photon index. 



The event rate of atmospheric z/-background is (see Anchordoqui et al. 2002b) 

= Aff / dE, Pu^^^iE,) Al^ioxio , (4) 



dN 
It 



where A^.^ is the effective area of the detector, Al^ioxio pa 3 x 10 sr is an assumed bin 
size for the observation (ICECUBE-generation detector, Karle et al. 2002), and d^B/dEy < 
0.2(E^/GeV)-3-2^ GeV-^ cm-^ s'^ sr"^ is the i/^ + atmospheric i/-flux (Volkova 1980, 
Lipari 1993). Here, P^^^{E^) 1.3 x 10^^ (£;^/TeV)°-^ denotes the probability that a i/ of 
energy E,y > 1 TeV, on a trajectory through the detector, produces a muon (Gaisser et al. 
1995). On the other hand, the i/-signal is 



dN 



Aff j dE,F,{E,)P,^^{E,) . (5) 



The signal-to-noise (S/N) ratio in the 1-20 TcV band is then {dN / dt\g) / dN / dt\^) = 
0.1 — 2.5 for one year of operation, depending on the photon spectral index. A NEMO- 
generation detector, with a bin size Afl ~ 0.3° x 0.3° (Riccobene et al. 2002), would have 
a S/N ratio roughly three times larger for the same integration time. The effect of neutrino 
oscillations may reduce the flux by a factor of three at Earth (Bilenky et al. 1999). However, 
for some possible values of the photon spectral index, this signal could be detected in a 
reasonable time span. The possible starburst neutrino spectra producing the previously 
quoted S/N ratio are plotted, together with an upper limit for the atmospheric neutrino 
background with vertical incidence, in Fig. Ib.^ 



average half the energy of the pion. Therefore, there is as much energy in photons as there is in neutrinos, 
and K = 1. 

^Regrettably, the south pole location of the ICECUBE telescope would render its observation of NGC 
253 nearly impossible; only neutrinos going through Earth can be identified. The smaller northern telescope 
ANTARES (0.1 km^, Becherini et al. 2002) would need much longer integration times to be able to distinguish 

NGC 253 from the background, if at all possible. Nevertheless, the analysis presented here can be directly 
applied to M82, a northern-sky starburst at the same distance as NGC 253. 
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6. Concluding remcirks 

Through hadronic interactions in their central regions, nearby starburst galaxies can 
generate high-energy 7-ray emission detectable at Earth. Proton illumination of the in- 
ner winds of massive stars in a high density CR environment may result in TeV sources 
without strong counterparts at MeV — GeV energies. The high number of early-type stars 
with strong stellar winds in this kind of galaxies might produce a strong collective effect, 
providing the setting for accelerating protons up to multi-TeV energies. In addition, we 
have found that starbursts might also be potential sources for km-scale neutrino telescopes 
within reasonable integration times, depending on the photon spectral index. The detection 
— or non-detection — of the neutrino signal may be an essential piece of evidence in the 
determination of the origin of the high-energy radiation. 

Gamma-rays from neutral pion decays will be produced not only in the stellar winds 
at the core of the starburst, but also in the diffuse matter among the stars. If we model 
the starburst region as a disk with a radius of ~ 300 pc and a thickness of ~ 70 pc, where 
the average density of H atoms in the diffuse ISM is ~ 300 cm~^ (Paglione et al. 1996), we 
get that the expected flux above 1 TeV is F{E > 1 TeV) ~ 3.9 x 10~^^ ph cm~^ s"\ for 
an average CR enhancement of ~ 1000 (see Torres et al. 2002 for details of calculation). 
This means an additional diffuse contribution of ~ 3 x 10^^ erg s~^ to the total TeV 7-ray 
luminosity of the starburst. In most of our models this diffuse contribution is only a small 
fraction of the total 7r°-decay 7-ray emission, which is dominated by interactions occurring 
in the stellar winds. 

Finally, it is important to remind that the CANGAROO observations are consistent with 
an extended source, whose size might be even bigger than the galaxy itself. If this is true, 
it would imply that there might be another component generating TeV radiation, related 
with the galactic halo. This component might probably be the result of Inverse Compton 
scattering of CMB photons by TeV electrons accelerated at the superwind. In such case, 
since the interactions occur in the Thomson regime, the spectrum should be harder than 
that produced by hadrons in the center of the starburst.® 

Comments by L. Anchordoqui, P. Benaglia, Y.M. Butt, C. Mauche, M.M. Kaufman- 
Bernado, M. Pohl, and an anonymous referee are gratefully acknowledged. G.E.R.'s research 
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ICET (under grant PIP N° 0430/98) and ANPCT (PICT 03-04881). The work of D.F.T. 
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We thank Prof. Yanagita for caUing our attention to this. 
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